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Abstract 

We study the flux noise in Josephson junction arrays in the critical regime 
above the Berezinskii-Kosterlitz-Thouless transition. In proximity coupled 
arrays a local ohmic damping for the phases is relevant, giving rise to anoma- 
lous vortex diffusion and a dynamic scaling of the flux noise in the critical 
region. It shows a crossover from white to 1 //-noise at a frequency lo^ cx £ _2 

with a dynamic exponent z = 2. 
PACS numbers: 74.40.+k, 74.50. +r 
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A variety of two dimensional systems undergo the Berezinskii-Kosterlitz-Thouless (BKT) 
pair-unbinding transition [|1],0. The transition between the high temperature (disordered) 
and the low temperature (coherent) phase is driven by the thermally excited vortices. These 
topological excitations form a two dimensional Coulomb gas. Below the BKT transition 
temperature Tbkt they are bound in vortex-antivortex pairs. Above the critical temperature 
the pairs dissociate and form a vortex plasma. In the plasma phase the vortex-vortex 
interaction is screened at a distance of the order of the correlation length £ which diverges 
at Tbkt- 

Arrays of Josephson junctions are prototype systems to study the BKT transition. Below 
the BKT transition the array is phase coherent and thus superconducting. On the contrary, 
in the plasma phase, free vortices destroy the coherence and the system is resistive, though 
each island remains superconducting. In the last decade there has been a great amount 
of work on the various aspects of the BKT transition in Josephson arrays (see Ref. [[3] for 
an overview). Experimental studies are based on electrical resistance Q, two-coil induc- 
tance and SQUID [7]|| measurements. 

Anomalous suppression of the vortex mobility has been observed in proximity coupled 
Josephson junction arrays by Theron et al. The vortex mobility vanishes logarithmically 
at low energies, showing that vortex motion in these arrays cannot be described by a Drude 
model with a constant mobility. This behavior has been explained either introducing a local 
ohmic damping for the phases of the superconducting wave function ||,[10| or by invoking 



long-range vortex interactions [11|. More recently a regime interpretation of the vortex 



dynamics was given in Ref. [12 



In a recent Letter, Shaw et al. [0] investigated the magnetic flux noise close to the BKT 
transition in an overdamped Josephson junction array by means of a SQUID. The flux-noise 
spectrum is defined as 

S*M = J dte^($(t)$(0)> , (1) 
where is the flux detected by the SQUID and is proportional to the vortex density 



integrated over the SQUID area. Near the transition the flux noise is white for frequencies 
uj < uj^ crossing over to a 1/ /-noise for to > u^. The crossover frequency vanishes as 
the transition is approached from above as uj^ oc , where £ is the correlation length and 
z ~ 2 is the dynamic exponent. All the data obtained at different temperatures collapse 
on a single scaling curve when plotted as a function of uo/uj^. The problem was studied by 



numerical simulations fl3|Jl4l| of both the Time Dependent Ginzburg-Landau (TDGL) and 
the Resistive Shunted Junction (RSJ) models. 

In this Letter we study theoretically the flux noise in the TDGL model near the BKT 
transition. We determine the scaling function of the flux noise analytically in the limit in 
which the SQID size is much smaller than the correlation lenght. 

In order to study the flux noise we map the XY-model with a local damping onto a 
2D (dynamic) Coulomb gas ||15|| . We determine the dynamic correlation functions using an 



improved Debye-Htickel approach for the Coulomb gas in the same spirit as in Ref. [[17 
The Euclidean action for the XY-model with local damping for the phases <p is 

m = * /; ir £ [1 - cosfeH - wW) ] + jfL dTdT ^ ( J^^ )' , (2) 

where Ej is the Josephson coupling between neighboring sites. The local ohmic damping 
introduces an interaction of the phase of a single grain at different times. In proximity 
coupled arrays, which consist of superconducting island on top of a metallic film, the local 
damping parameter a is related to the shunting resistance R of one island to the substrate, 
a = h/4e 2 R. 



A Villain transformation [18| is used to derive an effective dynamic action for the vortices. 



Here we follow the method described in Ref. [19] and obtain the dynamic action for the vortex 
degrees of freedom Viij) 
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The kernel D reads 



S[v] = \l drdr'^Vi^D^T - r') Vj {r') . (3) 

Z J0 i,j 



nil- \ An2Ej k 2 + 2\u„\/u a 
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with u a = 2nEj/a, and uj^ are the Matsubara frequencies. We choose the lattice spacing 
as the unit of length scales. In the static case this yields the standard 2D Coulomb gas. 

The evaluation of the flux noise is done by means of a two-step procedure in the spirit 
of Ref. ||17| . i)The physics taking place at length scales shorter than £ is taken into account 
by using the scaling properties of the flux noise, ii) At the scale of the order of the corre- 
lation lenght, all the dipoles have been integrated out and, since the physics above Tbkt is 
dominated by screening, a Debye-Hiickel approximation (with renormalized parameters) is 
used HO • 



The effect of bound pairs up to distances of the order of the correlation length can be 
account for by means of the scaling behaviour of the flux noise 

S 9 (u, I, T') = e z5 S^(ue zS , le~ s , T'(S)) . (5) 

At S* = ln(£/£o) all the vortices are integrated out up to a distance of the order of £. At 
this scale only free vortices are present and the Debye-Hiickel approximation can be used to 
calculate the r.h.s. of Eq. 

The effect of screening due to the presence of free vortices above Tbkt is analyzed in 
the Debye-Hiickel approximation. The Matsubara Green's function for the vortices in this 
approximation is 

G(k, = (vv) k ^ = Air2Ej pl D{kt(Jii) (6) 
The Kosterlitz correlation length £ diverges exponentially near the transition, £ = 



£o exp(&/y T' — T BKT ), with T' = T/Ej Pj7|p0|. An analytic continuation (1^1 — > —iu) 
yields the retarded vortex propagator G R {k,uj) = G(k,\u fl \ — > —iu) which is related to 
the spectral function of vortex density- density correlations by the fluctuation-dissipation 
theorem (vv)k )ti > — lm^G R (k,uj) 2T/u, for u <^T. This correlation function describes the 
equilibrium vortex fluctuations in the system and gives rise to the flux noise. 

Experiments on the magnetic flux noise detect the fluxes of the vortices from an effective 
area I 2 |6|-[8|]. The relevant quantity is the vortex density fluctuation integrated over the 



pick-up area $(t) = J t d 2 x v(x, t). After having performed the analytic continuation and by 
combining the scaling arguments (see Eq. (||)) with the Debye-Hiickel approximation (see 
Eq. (0)), we obtain the flux- noise spectrum 

s 9 (u) = *g /M^.^.W = ^^-f (* = ->v = nr)> ( 7 ) 

where the subscript DH means at at this scale the flux noise can be evaluated by means of 
the Debye-Hiickel approximation. The constant C takes into account the geometrical details 
of the experimental setup. We introduced the scaling function F 



F { X = ^ y = —)=- lm Jo dZ 



l + z- 
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(8) 



z — 2ix 

where characteristic frequency is L)£ = cg> (£o/£) 2 (we used a hard cutoff in the /c-space to 
integrate over the pick-up coil area). The dynamic exponent z = 2 directly emerges. The 
frequency scale u Q is related to microscopic parameters, u = 2irEj/(l; 2 a). The scale of the 
correlation length, £ 0) is of the order of the lattice spacing. Close to the transition (£ — > oo) 
the scaling function F reduces to 

, 1 + 3ZX , (\ — ix + Q(x) \ 7T . , 

F (x, y -> oo = Im— -— In - T 9 

v ' y 7 2Q(x) \l-ix-Q{x) J 4 v 7 

where Q(x) 2 = 1 — x 2 + 6ix. The flux noise in the relevant limits is given by 

( c$jj for ^ ^ 1 

W 5.H ~ 87r3 ?4 w < (io) 

The white noise occurs for frequencies u < u^, and 1/ /-noise for uo > uj^. The 1/ /-noise 
stems from the superposition of Lorentzian-shaped contributions at all length scales down 
to the dimension of the SQUID. For finite £ the 1// noise crosses over to a l// 2 -noise. The 
l//-noise is intermediate in the frequency range 1 < uj/uj^ < AttC, 2 /1 2 . Close to the transition 
this range exceeds several orders of magnitude. 

We finally apply our results to the experiment by Shaw et al. |7| . The critical coupling in 
the experiment is Tg KT = 0.06, thus the bare fugacity at criticality is yo = 0.35. Assuming 
that the bare fugacity is only weakly temperature dependent, the numerical iteration of 



the RGE gives to a good approximation T'(S) oc exp(45), as upon renormalization the 
limit T'(S) ^> 2/tt is reached, where y(S) oc exp(25) (this scaling behavior reflects the 
fact that the vortex-pair density is constant). Thus T'(5*)/£ 4 is temperature independent 
and we can identify our scaling function F with the experimental scaling behavior. These 
last considerations should be considered more as additional assumptions whose validity is 
mostly based in the comparison with the experiments which can be made quantitative. The 
characteristic frequency was found to be fo = 2.1 x 10 6 Hz. With the estimate Ej = 10K 
we can extract = 10 6 , corresponding to ohmic shunts to the ground with a resistance 
R = 6 mf2 x £q . We expect £o to be of the order of one lattice spacing. Yet another possibility 
is to determine the resistance R and deduce the value of £ - The effective shunting resistance 
of one island to infinity can be estimated be of the order of 10 mfl. This also yields the scale 
of the correlation length to be of the order of one lattice spacing. 

In summary we calculated the flux noise above the BKT transition and derived its scaling 
behavior analytically in the case where the correlation length exceeds the size of the pick-up 
loop. We find a crossover from white to l//-noise. We related the crossover frequency to 
microscopic parameters of the array. Our results may explain a recent experiments by Shaw 
et al. 0. 

We are indebted to T.J. Shaw for useful conversations and for helping us in comparing 
our work with the experiment. We would like to thank J. Clarke, D.-H. Lee, P. Minnhagen, 
Gerd Schon, A. Schmid, and M. Tinkham for fruitful discussions. This work was supported 
by "SFB 195" of the "Deutsche Forschungsgemeinschaft" and the EU-grant TMR-FMRX- 
CT-97-0143. 



6 



REFERENCES 

[1] V.L. Berezinskii, Zh. Eksp. Teor. Fiz. 59, 207 (1970) [Sov. Phys. JETP 32, 493 (1971)]. 

[2] J.M. Kosterlitz and D.J. Thouless, J. Phys. C 6, 1181 (1973); J.M. Kosterlitz, J. Phys. 
C 7, 1046 (1974). 

[3] For a number of references see: Proceedings of the NATO ARW on Coherence in 
Superconducting Networks, edited by J.E. Mooij and G. Schon, Physica B 152,1-302 
(1988); Proceedings of the ICTP Workshop on Josephson Junction Arrays, edited by 
H.A. Cerdeira and S.R. Shenoy, Physica B 222, 253-406 (1996). 

[4] D.J. Resnik, J.C. Garland, J.T. Boyd, S. Shoemaker, and R.S. Newrock, Phys. Rev. 
Lett. 47, 1542 (1981); D.W. Abraham, C.J. Lobb, M. Tinkham, and T.M. Klapwijk, 
Phys. Rev. B 26, 5268 (1982); H.S.J, van der Zant, H.A. Rijken, and J.E. Mooij, J. Low 
Temp. Phys. 82, 67 (1991). 

[5] Ch. Lehmann, Ph. Lerch, G.-A. Racine, and P. Martinoli, Phys. Rev. Lett. 56, 1291 
(1986). 

[6] R. Theron, J.-B. Simond, Ch. Lehmann, H. Beck, P. Martinoli, and P. Minnhagen, 
Phys. Rev. Lett. 71, 1246 (1993). 

[7] T.J. Shaw, M.J. Ferrari, L.L. Sohn, D.-H. Lee, M. Tinkham, and J. Clarke, Phys. Rev. 
Lett. 76, 2551 (1996). 

[8] Ph. Lerch, Ch. Lehmann, R. Theron, and P. Martinoli, Helv. Phys. Acta 65, 389 (1992). 
[9] H. Beck, Phys. Rev. B 49, 6153 (1994). 
[10] S.E. Korshunov, Phys. Rev. B 50, 13616 (1994). 

[11] J. Houlrik, A. Jonsson, and P. Minnhagen, Phys. Rev. B 50, 3953 (1994); K. Holmlund 
and P. Minnhagen, Phys. Rev. B 54, 523 (1996). 

[12] M. Capezzali, H. Beck and S.R. Shenoy, Phys. Rev. Lett. 78, 523 (1997). 

7 



[13] P.H.E. Tiesinga, T.J. Hagenaars, J.E. van Himbergen, and J.V. Jose, Phys. Rev. Lett. 
78, 519 (1997). 



[14] I.-J. Hwang, S. Ryu and D. Stroud, |cond-mat 97U41U8 



[15] P. Minnhagen, Rev. Mod. Phys. 59, 1001 (1987). 

[16] B.I. Halperin and D.R. Nelson, J. Low Temp. Phys. 36, 599 (1979). 

[17] This procedure was developed to study the specific heat and the superfluid density 
above T BK t in A.N. Berker and D.R. Nelson, Phys. Rev. B 19, 2488 (1979). 

[18] J. Villain, J. de Physique 36, 581 (1975). 

[19] R. Fazio and G. Sch6n, Phys. Rev. B 43, 5307 (1991). 

[20] This exponential temperature dependence of the correlation length may, in a strict 
renormalization group sense, not directly reflect the critical BKT-behavior, as pointed 
out by P. Minnhagen and P. Olsson, Phys. Rev. B 45, 10557 (1992). 



8 



